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We explore the phenomenological quark-lepton mass matrices, which are devised fol- 



lowing the S3 flavour symmetry principle, yet fully consistent with SU(5) gauge models 
with the Higgs particles of 5, 45 and their conjugates. The model contains 10 free param- 
eters altogether. When 6 parameters are fixed by charge 2/3 quark and charged lepton 
masses, charge(— 1/3) quark masses and all quark mixing matrix elements are predicted to 
be close to experiment, leaving some narrow ranges still as freedom. Further specification 
of 3 more parameters (using, m^, Cabibbo angle and \ V c b\) suffices to fix the quark mixing 
matrix nearly completely, and all elements come out to be in accurate agreement with 
experiment. We obtain the CP violation phase as a prediction. 



Grand unification of gauge theories has not given much hint to our understanding of 
the quark-lepton mass spectrum. The most successful among many attempts is perhaps 
the "prediction" of the mass relations between the charge —1/3 quarks (referred simply to 
as down quarks) and the charged leptons by an introduction of 45-plet Higgs in addition 
to the standard 5-plet within the SU(5) grand unification IJ. The relation reads, 

m d = 3m e , m s = ^m M , m b = m T , (1) 

which is often called Georgi-Jarlskog mass relation. No successful prediction, however, 
has been known for the charge 2/3 quark spectrum and hence for quark mixing. Even 
more difficult is to understand the CP violation phase and its origin. 

The only approach which turns out to be "successful" in giving correct mass-mixing 
relations for quarks is an empirical approach, where some discrete symmetry is imposed on 
the form of mass matrices and fix parameters using some quark masses as input & 0, ffl . 
One of the most unsatisfactory aspect of such approaches was that its consistency is not 
clear with the unified gauge model, which anyway we must impose at some level. If one 
would impose the compatibility with a gauge model in a straightforward way, we are 
usually led to unwanted relations for quark and lepton masses as a remnant of prototype 
unified gauge models. 

We have devised in a previous paper || phenomenological quark-lepton mass matrices 
based on the S3 permutation symmetry principle in a manner fully compatible with SU(5) 
grand unification. This model results in an approximate Georgi-Jarlskog relation and a 
mixing angle pattern for the quarks, which are in decent agreement with experiment. 
The model also successfully applies to the neutrino mass-mixing problem with bimaximal 
mixing as a natural outcome ||. For sake of simplicity of the argument, we have assumed 
in || the matrix elements to be all real, ignoring all phases which could in principle appear 
therein; also, accurate agreement with experiment was not sought for the mixing angles. 
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In this paper we extend our analysis allowing for full degrees of freedom of the mass 
matrices, trying to cure the defects of the previous model. We have now more parameters, 
but they are still tightly constrained within the model, and this leads to a prediction of 
the CP violating phase. Here we should quote an earlier work of Fritzsch [[!]], who has 
also derived the CP violating phase in his matrix model approach. 

We begin with the Yukawa coupling in the SU(5) model: 



^Yukawa = Y{h H )uij^i^H + ^(45^10*1(^45^ 

+Y(5* H ) D/E i j 5*10 j 5* H + F(45^ /Eii 5* 10 i 45 tf + K{h H h H ) vij h*5\ 
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(2) 



where bald face symbols with suffix H denote Higgs scalars of a specified multiplet, and 
those with suffix i or j (refer to flavour) are SU(5) matter fields, 5* = (d{, d\, d\, e~ , v e )Li 
and 10j = (ul, ...,Ui, •••,c?f, ...,e + )ij- We specify down quarks and charged leptons with 
suffix D/E as they are unified. The last term of eq.(|2|) is an effective neutrino coupling 
where the neutrino is assumed to be of the Majorana type. We suppose that it is induced 
from heavy Majorana right-handed neutrinos N{ (SU(5) singlet) with mass Mr ||, so 
that 45#45# does not appear in (fj). 

We postulate the mass matrices of the form 0: 
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for the down-quarks, and 
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for the charged leptons. Here, the main part of the mass matrices is induced by a 5 plet, 
which is S3 permutation symmetry invariant. We write it So where 10 and 5 

refer to representations of fermions. We break S3 symmetry in a hierarchical manner. We 
introduce 5 terms to break S^xSg down to S^xSf . We assume that 5d and 5e elements 



are generated from the coupling to a 45*#-plet Higgs scalar. This assignment removes 
the unwanted down-quark charged lepton mass degeneracy of the minimal SU(5) model, 
but produces the Georgi-Jarlskog mass relation. Further symmetry breaking is caused 
by e terms (e <C 5) in a way consistent with 5*h to allow further adjustment of mass 
hierarchy. We have then 

e E = e D , 5 E = —35 D . (5) 

For the up quark masses, respecting the same principle of the symmetry breaking pattern 
as with the down quark/charged lepton sector), we write 
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(6) 



where the main term and e terms arise from 5# and 5 terms from 45 h- Note that the 
texture of 5u in eq.(^) is an unique invariant of S^xSf among the anti-symmetric (45) 
mass matrix. 

In general, parameters €d, #d, £u and 5jj are complex, and we can express them as 
Q = \ti\e iai and 5i = \6i\e lf3l (i = U, D), whereas we assumed them to be all real in EJ. We 
write the matrices in the hierarchical base M by applying an orthogonal transformation 
in order to investigate the structure of the phase, 
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(9) 
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Extra phases in My and My are removed by a phase transformation applied to the 
quark fields qy and q^. With diagonal phase matrices Py and Py, we obtain 
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P D M D P D ~ ^ 
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and 
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where 
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In eq. fllPp , only leading terms are retained in the (2,2), (2,3), (3,2) and (3,3) elements 
for simplicity of our expressions, while we keep all terms when we carry out a numerical 
analysis. The expression of eq.(|TT|) is exact. Phase matrices Py and Py contribute to the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix as: 



VcKM = UyPyP D U D , 



(13) 



where Uy and Uy are unitary matrices, which diagonalize My and My. For convenience, 
we define the phase matrix Q by 



Q 



PuPd 



where a = (ay — Py) — (ay — (3y) and 
The unitary matrix Uy is given by 
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where s? = sin 6P and cP = cos$P with 
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The down-quark masses are then, 
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(17) 



The charged- lept on masses are given as 



K D {1 - -\8 D \ cosf3 D ), 
--K D \S D \ (l + -\5 D \ cos/3 D j , 

in,,,, ~ -K DTT - r ( 1 — -t|0d| cos pd 1 ■ 
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The matrix for the up quark sector Uu is obtained as 
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(19) 



where 
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2 \5 V \ 
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The up- quark masses are 




c s 17 

> A 3 
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(20) 



m t ~ ^ ( 1- 



--Kjj- 



(21) 



6 u |^| 2 

Substituting [//, and Ujj into eq . (|i~3"f ) , we obtain the CKM matrix in terms of the 



quark masses and the phase parameteres. After taking c\ ~ c % 2 — c| ~ l(i = D, U) and 
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neglecting s\s\, s l 2 s\ and s\s\s % 3 (i = D,U) terms, we obtain the CKM matrix elements 
approximately as 



v ud * Ji - 

V m x 



V us * J^ e *(^M - pie* 
V rn s V m c 



V ub ~ _ /_i_V(*-A>) - —=—J—e %T + J— + W 17 

2v2y \/2m b \ m c ]] m t y 

V m.= V m fi 
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V2 m fe y ra t 



yr J_ / m d m s | 1 m s \ m d m c ^ T l 2m u m c ci(a _ 0D+0u) 

2v / 2y m 6 y/2m b y m s y m s m t y mf 



1 m s /m^ iT 



V2 m 6 y m t 

V tb ~ 1. (22) 
For the CP violation parameter, the angle 7 that enters in the unitarity triangle || reads, 

j~-(a-(3 D )+sm (-—, = - nrr^ -sin I J — - , 23 

V 2V2 m b \V 13 \ J Vy m c \V 12 \ J 

in the case of a small r. 

These expressions agree with the result of a numerical calculation (without any ap- 
proximations) within 10% error. Although there are four phase parameters a, r, (3d and 
Pu in our matrices, the CKM matrix is determined in practice by only two of them, a — (3d 
and r, because the last terms of V ub and V td , which contain phases other than the two, are 
strongly suppressed compared with other terms. In |J we have shown that these matrix 
elements, where all phases are completely dropped, yield a resonable description for all 



CKM matrix elements. The point of the present paper is that all CKM matrix elements, 
including the CP violation phase, are completely determined, once the above two phases 
are fixed by the adjustment of two of the CKM matrix element, say V us and V c b- 

While proceeding to a numerical analysis, we need some care as to the input data. 



Since the CKM matrix elements in eq . (|22|) and masses eqs.f pOjpTD are discussed at the 
SU(5) GUT scale, predictions should also be compared at the GUT scale, rather than 
at the electroweak scale. Since a supersymmetric extension is the only way to make the 
SU(5) GUT viable, we take quark and lepton masses at the GUT scale obtained in the 
minimal SUSY model (MSSM) with the aid of renormalisation group equation (RGE) 
incorporating two-loop|[fO||. These mass parameters are given in Table 1. 



Let us first discuss charge —1/3 quark and charged lepton masses. With the three 
charged lepton masses as input, the parameters in the charged lepton/down quark sector 
are determined to be Kd = 1.203GeV, \5d\ = 0.080 and \eo\ = 0.0104, which in fact 
satisfies K D ^> \5d\ ^> \^d\- The ratio of <i-quark mass to electron mass is given, using 
eqs.(0) and flU), by 

m d . o 1 + j\$d\cos(3 d ^ 





cos (3d 


i-I 


S D 


cos (3d 



This shows the dependence of |m<i/m e | on further two parameters. By carrying out an 
accurate numerical calculation with \5d\ = 0.08, we see that this ratio takes a value 
between 2.3 and 3.7, which corresponds to = 0.76 — 1.20MeV, when (3d is varied form 
to +7T. This is compared with the "experimental value", = 1.3 ± 0.2MeV (i.e., the 
ratio is 3.4-4.6). Requiring an agreement with experiment leads to 0.7 < cos (3d- Here we 
take cos (3d = 1 for further analysis. For this case we have rrid = 1.20MeV. The prediction 
for other down quark masses is given in Table 1. The values of rrib and m s are somewhat 
larger, but taken as acceptable when we consider the uncertainty in the mass analysis 
using low-order perturbation theory. 

Now we are concerend with the CKM matrix element. We obtain the three parameters 
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of the up quark sector to be K u = 129GeV, \5u\ = 0.1026 and \e v \ = 0.00071 (K u > 
\$u\ ^> \ e u\ being satisfied) if the central values in Table 1 are adopted for up-quark 
masses. We note, however, a large uncertainty in the estiamte of the top quark mass at 
the GUT scale, arising from the fact that the top quark mass is near the fixed point of 
RGE at the electroweak scale. The prediction of the CKM matrix depends on top quark 
mass used as input; we take account of this large error, 89 — 325GeV (for a given pole mass 
mt (pole) = 180 ± 12GeV), for further analysis of the matrices. As with the case of the 
mass, we must take account of the running effect of the CKM matrix elements. It is known 
that the elements V u d, V us , V c d, V cs and V t b are nearly constant during running between the 
electroweak and the GUT scales. On the other hand, all others are affected by the large 
Yukawa coupling of the top quark by 10-20%. We may use |K S | = 0.217 - 0.224 (GUT 
scale) to constrain the phase = o — r — [3d- We then obtain = (±60.5°) — (±68.5°). 
Another phase r is fixed by |V^,| = 0.0347 — 0.036; we obtain r = — (±22)° for m t = 
129GeV at the GUT scale. All parameters are now fixed, and a prediction is given for the 
full CKM matrix elements, as presented in Table 2. In this table we also give experimental 



values and estimates at the GUT scale after running [0. All predictions are within the 
uncertainty of the empirical values. The agreement with experiment is not spoiled even 
if we take the upper limit value m t = 325 GeV at the GUT scale. On the other hand, if 
we decrease the top quark mass the agreement is disturbed for \V c b\: the prediction goes 
out of the upper limit of the range allowed by experiment, 0.030 — 0.036. This limits our 
consideration to 123GeV <m t < 325GeV. 

The most interesting prediction is perhaps that of the CP violation parameter. We 
obtain 

7 = 106° - 114°, (25) 
for m t = 129 GeV. This is equivalent to the vertex position of p = —0.098 ~ —0.086 
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and 7] = 0.221 ~ 0.309 in Wolfenstein's p — r\ plane [II]. [] If the range of the input 
m t is relaxed as discussed above, the prediction for 7 becomes 76° — 114°. This value 
is similar to the prediction of Fritzsch 0, 7 = 72° ~ 76°. In Table 3, we present the 
prediction of the several key parameters for various top-quark masses. Here, we take 
iKtfc/'Ve&l as one of the indicators. Fig. 1 is a summary of our prediction for the CKM 
matrix in the p — rj plane, where currently available experimental constraints are also 
plotted: (i) e K parameter with B K = 0.6 - 0.9 (ii) \V ub /V cb \ = 0.08 ±0.02 and (iii) AM Bd 



with ^JB Bd F Bd = 200 ± 40MeV. 

We have explored the full content of the phenomenological quark-lepton mass matrices, 
which are derived from S3 flavor symmetry principle and its hierarchical breaking in the 
framework of SU(5) gauge models, proposed in |J. Our mass matrices have 10 free 
parameters altogether, 6 real K iy \5i\, \ei\(i = D, U) and 4 phases (3d, &, Pd, t, excluding 
neutrino parameters. Using (m u , m c , m t ) and (m e , m M , m T ) as input, we have a prediction 
in decent agreement with experiment for the rest of the physical parameters, irrespective 
of the free parameters left unspecified. Requiring more precise agreement for irtd and 
two mixing angles, say \V US \ and \V c b\, we have a complete determination of the CKM 
matrix elements including phases. The obtained matrix shows an excellent agreement 
with experiment within the present experimental accuracy. The most interesting among 
others is the prediction of the CP violation phase, which would soon be tested to higher 
accuracy with B-factories. 

We have also studied the mixing problem for the lepton sector including neutrinos in 
a way parallel to the quark sector, but have found little to add to the previous work of 
ref.0, except for phases. With the allowed parameter range the net CP violating phase, as 

defined by Jarlskog fL2"| , is as small as Jqp < 10~ 4 because the CP violating phase enters 

lr These predictions should in principle be compared with the corresponding values at the GUT scale. 
Fortunately, the energy scale dependence of these quantities is very weak |l(J , and we can safely neglect 
the running effect. 
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only the symmetry breaking terms. This is too small to arouse any phenomenlogical 
interests. 
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m e 
(MeV) 


(MeV) 


m T 
(GeV) 


m d 
(MeV) 


m s 
(MeV) 


m b 
(GeV) 


m u 
(MeV) 


m c 

(GeV) 


m t 
(GeV) 


"exp." 


0.325 


68.60 


1.171 


1.3 ±0.2 


26.511;$ 


1.00 ±0.04 


1.0 ±0.2 


3021P 


129l^ 6 


pred. 


input 


input 


input 


1.20 
(input) 


19.6 


1.22 


input 


input 


input 



Table 1: Input quark- lepton mass parameters and the prediction of our model at the GUT 
scale, rrid is the value with cos (3d = 1- 





"exp." 


value at GUT scale 


prediction at GUT scale 


\v ud \ 


0.975 - 0.976 


0.975 - 0.976 


0.975 - 0.976 


\v us \ 


0.217-0.224 


0.217-0.224 


0.219-0.221 


\v ub \ 


0.0018 - 0.0045 


0.0015- 0.0040 


0.0019- 0.0026 


\Vcd\ 


0.217-0.224 


0.217-0.224 


0.219 - 0.220 


\Vcs\ 


0.974 - 0.976 


0.974 - 0.976 


0.975 


\Vcb\ 


0.036 - 0.042 


0.030 - 0.036 


0.035 - 0.036 


\Vu\ 


0.004- 0.013 


0.0035 -0.011 


0.007- 0.008 


\Vu\ 


0.035 - 0.042 


0.030 - 0.036 


0.034- 0.035 


\V tb \ 


0.999 


0.999 - 1.000 


0.999 



Table 2: The CKM matrix elements. The first column shows experiment, the second is 
the values estimated at the GUT scale. The third column is the prediction of our model 
with m t = 129GeV at GUT scale. The underlined values are input. 



m t (GeV) 


M(°) 


\v cb \ 


\v uh /v cb \ 


P 


V 


7(°) 


129 


0-22 


0.035-0.036 


0.055-0.072 


-0.098 - -0.086 


0.221-0.309 


106-114 


150 


0-43 


0.031-0.036 


0.056-0.085 


-0.116 - -0.045 


0.216-0.377 


97-118 


200 


41-69 


0.030-0.036 


0.091-0.094 


-0.068 - -0.033 


0.402-0.421 


86-100 


250 


60-85 


0.030-0.036 


0.096-0.099 


0.004 - 0.082 


0.425-0.445 


79-90 


325 


76-90 


0.030-0.033 


0.099-0.103 


0.068 - 0.105 


0.435-0.457 


76-82 



Table 3: Predictions of the CKM matrix elements and the unitarity triangle. m t is the 
value at the GUT scale. 
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Figure 1: Predicted vertex of the unitarity triangle on the p — r] plane for m t (GUT) 
varying between 123—325 GeV. Experimental constraints from e K) \V ub /V cb \ and AM Bd 
are also plotted. 
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